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Based on the degeneracy of the d zx and d yz orbitals in 
S^RuCU it is argued that the Cooper pairs condense in or- 
bital singlets. Together with the spin-triplet wave functions 
the real-space wave function then is symmetric. Considering 
interaction effects the order parameter is found to have A\ g 
symmetry consistent with a number of experimental observa- 
tions. The sensitivity of the material on non-magnetic impu- 
rities follows in a straightforward manner from the orbital- 
singlet configuration. 



With the discovery of the high temperature supercon- 
ductors a whole class of transition metal oxides became a 
focal point in condensed matter research. These materi- 
als exhibit many unconventional properties whose inter- 
pretation has so far generally proved controversial. An 
example that attracted a lot of attention is Sr2Ru04. Rs 
normal state properties are Fermi liquid like in the tem- 
perature range T c < T < 30 K |[§ but below T c < 1.5 
K the material is an unconventional superconductor || 
since a number of experimental probes J^— 16| show that 
the paired electrons carry a magnetic moment. In spite of 
the large interest that the superconductivity in Sr2Ru04 
has attracted an unambiguous understanding of the elec- 
tronic correlations has not yet evolved. 

Rice and Sigrist |?J proposed that the superconducting 
order parameter has p-wave symmetry promoted by fer- 
romagnetic correlations by analogy with 3 He. This idea 
is supported by experiments that show that the static 
magnetic properties of Sr2Ru04 are the same in the nor- 

However, there 
for the p-wave 



mal and the superconducting phase [ 4 
is no conclusive experimental proof 
symmetry of the superconducting order parameter and 
no indications of ferromagnetic correlations have been 
found either in neutron scattering investigations [ fll| or 
other approaches fl^-[T^|- 

Furthermore, the specific heat fi~5|| , nuclear quadrupole 
resonance (NQR) ||, and thermal conductivity Jll| are 
consistent with two-dimensional gapless fluctuations in 
the superconducting phase of Sr2Ru04, which are incom- 
patible with the analogy to superfluid 3 He. One possible 
scenario is the existence of line nodes similar to those 
in the superconducting cuprates ]l~7p . Since vertical line 
nodes have been ruled out by thermal conductivity mea- 
surements horizontal line nodes in the subsystem 
of the d zx and d yz electrons have been proposed Jl8| . 
One weakness of the latter picture is that it requires the 
fine tuning of various interaction strengths |19[| , while no 



double gap structures have been observed in Andreev re- 
flection spectroscopy data ||. 

In this letter it is shown how the degeneracy of the 
Ru 4+ d zx and d yz orbitals allows for a straightforward 
description of the unconventional superconductivity in 
Sr2Ru04 that is consistent with the experimental obser- 
vations. The possibility of mixed orbital pairing leading 
to S = 1 spin-triplet Cooper pairs through Hund's rule 
coupling has been raised implicitly by Baskaran |2(]] . The 
"active" d zx and d yz orbitals drive the superconducting 
instability because they have the larger inter-plane elec- 
tronic overlap |2l| l . This is supported by the recently im- 
plied increase of T c upon uniaxial pressure [2^| along the 
crystallographic c axis since the inter-plane coupling is 
increased. Such pairing is umklapp scattering enhanced 
by the body centered tetragonal lattice |p3| . 

The results of the approach can be summarized as fol- 
lows. The Cooper pairs form orbital singlets allowing 
for an even parity real-space wavefunction in spite of the 
spin-triplet configuration. Taking into account the rel- 
atively strong interaction effects in the system [|| this 
allows for an almost homogeneous gap function, consis- 
tent with the experimental observations P-flO 24 1. Any 
impurity or defect [ p5| , p6[ locally breaks the symmetry of 
the d zx and d yz orbitals and thus acts as a pair breaker 
in strict analogy to magnetic impurities in a spin-singlet 
superconductor pT| |. The quadratic temperature depen- 
dence of the specific heat [[Hi] follows from fluctuations 
of the internal degrees of freedom of the order parame- 
ter p3| ]. On the other hand, the pair correlations for the 
d xy electrons are induced by the interband proximity ef- 
fect. Since this effect is usually strong Jl^] a single gap is 
assumed leading to consistency with Andreev reflection 
experiments (§). 

In the subspace of the degenerate d zx and d yz orbitals 
the possible order parameters can be classified in stan- 
dard notation |28| as orbital-singlet spin-triplet compo- 
nents 



t\ 



t,fJ, 



E 



4 c 1 * , ,) 



and orbital-triplet spin-singlet components 

/pt,rf\ = a v ct m a v 



(i) 



(2) 



Here fi = x,y,z labels the triplet components, are 
Pauli matrices, and v a are the usual electron creation 
operators on site n in orbital d zx (y — x) or d yz (y = y) 
with spin a. 
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FIG. 1. Product of the an- 

gular components d zx (r /r)d yz (r /r) — d zx (k/k)d yz (k / k) in 
the real or Fourier space projection of the pair wave function 
Eq. §. 



In the presence of Hund's rule coupling the spin- 
triplet states are energetically favored over the spin- 
singlet states and thus form the ground state of the 
Cooper pair condensate. Possible effects from spin-orbit 
coupling A w 0.1 eV Q are over-compensated by the 
larger Hund's rule coupling Jh ~ 0.2 — 0.4 cV Q since 
A<J H @. 

Following Eq. (Q) the condensate wavefunction \ip^) = 
P^^jO) can be factorized into spatial, spin-triplet, and 
orbital-singlet y) s ) contributions. Since the orbital 
singlet is odd and the spin triplet is even under electron 
exchange the real space projection of the wave function 
has to be even. This requirement is satisfied for electrons 
in the d zx and d yz orbitals, for which the total wavefunc- 
tion is Bl: 



(R n + r\tp s ) = d zx (r) d yz (r) \x, y) s | t, I), 



(3) 



R n are the real-space coordinates of the n th Ru ion, r 
are the coordinates relative to Rn. The product of the 
angular components d zx (r /r)d yz (r /r) has even parity as 
shown in Fig. |l|. 

The most striking evidence for orbital-singlet pairing 
in Sr2Ru04 is the sensitivity of the superconductivity to 
impurities [^5| and crystal defects |2(J , which can be un- 
derstood by analogy to the effect of magnetic impurities 
in a spin-singlet superconductor. In the latter the mag- 
netic impurities locally break spin-rotational invariance 
and thus act as pair breakers for spin-singlet Cooper pairs 
|p7f . Similarly, any impurity — magnetic, non-magnetic, 
or crystal defect — locally breaks the rotational symme- 
try of the lattice and thus the symmetry between the 
d zx and d yz orbitals. Consequently impurities are pair 
breaking in the orbital-singlet superconductor described 
by Eq. ([}]) in strict analogy to magnetic impurities in a 
spin-singlet superconductor. The resulting quantitative 
applicability of the theory of Abrikosov and Gor'kov to 
Si'2Ru04 is impressively demonstrated in Refs. 25 2t|. 

As a next step it is necessary to study the supercon- 
ducting gap function in order to interpret the numerous 
directionally dependent experimental probes. The gap 
function = (Afc iX , Afc iJ( , Afc j2 ) is given in the Fourier 



representation of Eq. ([!]) with (P'l^) = J2k ^fc vm 



(7,£T ,<T" VM' 



(4) 

It is determined in principle by solving the Eliashberg 
equations 



T 



Vc 



0.2 



j Vfa] e~ SsG ^ (5) 



self-consistently. Vb,2 is the effective pairing potential. 
The action Ssg has been derived using the quasi one- 
dimensionality of the kinetic energy of the d zx and d yz 
electrons and includes the intermediate coupling on-site 
interaction non-perturbatively [32] |. It depends on the 
four Bose fields, which can be considered as charge, fla- 
vor, spin and spin-flavor fields in analogy to the two- 
channel Kondo problem Q , and includes mass generat- 
ing terms in the superconducting state p3| ]. The treat- 
ment of such a four-component, two-dimensional sine- 
Gordon action is quite involved and is only possible using 
approximations. 

However, to investigate the gap function in Sr2Ru04 
it is sufficient to apply qualitative physical arguments. 
Starting with the investigation of the wavefunction sym- 
metry within the non-interacting, local picture and then 
analysing the expected influence of strong interactions it 
turns out that a rather homogenous gap function must 
be expected. 

To establish the wavefunction symmetry in momentum 
space it is useful to write [p4[ 



exp (ikr) — 4ir 



Fdkr) 
kr 



YC m {k/k)i l Y lm {r/r), (6) 



so that the angular components in real and Fourier space 
factorize. Fi(kr) is a regular spherical Bessel function 
and does not depend on the magnetization quantum 
number m. Since the d zx and d yz orbitals are linear 
combinations of the orthogonal spherical harmonics Y 2 ±i 
the angular part of the pair wavefunction projection onto 
Fourier space, (klip'p, has the same symmetry as in 
real space, i.e., d zx (r/r)d yz (r/r) = d zx (k/k)d yz (k/k) 



[Fig. | 

k k 

y ' y 



Introducing the rotation operator 71^/2 '■ k x 
-> — k T one has 



TZ v/2 d zx (k/k)d yz (k/k) = -d zx {k/k)d yz {k/k) , (7) 
K 7T / 2 \x,y) s = -\x,y) s , (8) 

and U„/ 2 \ T, Dm = I T,4>m- Consequently 

ft*/ 2 <fcK) - (*#*)■ (9) 

In other words the wavefunction is even under a rotation 
of 90° since both angular and orbital-singlet contributions 
are odd under that rotation. We therefore expect the gap 
function to be of extended s-wave symmetry. 
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FIG. 2. Representation of the Fermi surfaces formed 
by the d zx and d yz bands j3(jy. Dashed lines: idealized 
one-dimensional band. Full lines: hybridized bands. Dots: 
points at the Fermi surface where electrons can pair with op- 
posite momentum. Small black dots: idealized ID, larger 
dots: hybridized. 



In a group theoretical context the six possible pairing 
states described by the pair operators in Eqs. ([!]) and 
(0^) find their analogies in the possible pairing states of 
the tetragonal point group D^h J35|. Since the angu- 
lar Fourier space part d zx (k/k)d yz (k/k) has even par- 
ity and the pair wavefunction is invariant under rota- 
tion of 90° the state with either A\ g or A 2g symme- 
try must be realized. It is usefull to define the mir- 
ror operators A4 X : y — > —y and M y : x — > —x 
as well as M. x ■ V —> —y and M y : x — > —x with 
x = (x + y)/ \/2 and y = (x — y)/ \f2. Note that x and y 
define a reference frame rotated by 7r/4. Applying these 
to d zx (k/k)d yz (k/k) and \x,y) s reveals the An g symme- 
try of 

However, in the real system the symmetry of the gap 
function will be significantly altered by hybridization 
and — more importantly — interaction effects. Starting 
out by considering the non-interacting case the orbital- 
singlet superconducting instability can be formulated fol- 
lowing BCS p3| and electrons with opposite momentum 
can only pair at the four points (small black dots in 
Fig. ||) of the Brillouin zone where the idealized, one- 



dimensional Fermi surfaces 
cross (dashed lines in Fig. E 
the d zx and d 



!6l of the d zx and d yz bands 
In a more realistic picture 
and 



yz bands are weakly hybridized |3q 
form the a and (3 sheets of the Fermi surface (full lines in 
Fig. ||) . It can be shown that then only the eight points 
on the a and f3 sheets indicated by the larger dots in Fig. 
^ contribute to the pair formation. Such a small phase 
space for the pairing is consistent with the A 2g symme- 
try discussed above but would be inconsistent with the 
large specific heat anomaly at the superconducting phase 
transition M |. 

That the neglect of the interaction clearly represents 
an unjustified oversimplification of Eq. (|^) becomes ob- 
vious from the significance of the on-site interactions ]32] ] 
for the observed |]ll| , |3"7f strong magnetic in-plane corre- 
lations. An estimate of how the interactions increase the 
pairing phase space is possible by noting that the dom- 
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FIG. 3. Representation of the idealized one-dimensional 
Fermi surfaces formed by the d zx and d yz bands (dashed 
lines), (a) Magnetic momentum transfer (2&f, 2&f) and three 
combinations with reciprocal lattice vectors, (b) and (c) show 
the resulting momentum transfer (double arrows) that allows 
for mixed orbital pairing on many points of the Fermi sur- 
faces. The corresponding momenta of the electrons forming 
a pair are illustrated by the black arrows starting from the 
zone center in panel (b). (d) shows the momentum transfer 
for (4fc F , 4fc F ). 



inant magnetic correlations can be described as gapless, 
quasi one-dimensional fluctuations at momentum trans- 
fer q i = (±2/cF,±2fcF) modulus a reciprocal lattice vec- 
tor | p2[ . The arrows in Fig. |3|(a) show the momentum 
transfer of q 1 = (2/cF,2fcp) and three combinations with 
reciprocal lattice vectors Go = (0, 2ir), G\ — (2ir, 2tt), 
and G2 = (2tt, 0) in units of the reciprocal lattice spac- 
ing 1/a. 

The back-scattering terms in the action of E q. (H ) cou- 
ple magnetic and charge degrees of freedom [ p2[ 7 The 
Cooper pairs can thus scatter elastically off the gap- 
less magnetic excitations modulus any reciprocal lat- 
tice vector Gi, i.e., {c-k,y,<T'Ck±q i ±G j ,x,(r) ^ as indi- 
cated by the black arrows in Fig. ||(b). The resulting 
momentum transfer allows for mixed orbital pairing on 
many points of the Fermi surfaces formed by the ideal- 
ized one-dimensional d zx and d yz bands as indicated in 
Fig. ||(b) and (c). Including also higher order contribu- 
tions allows for an even more homogeneous distribution 
of paired electrons across the Fermi surfaces as indicated 
for (±4fc F ,±4fc F ) in Fig|(d). 

This qualitative discussion shows that interactions can 
be held accountable for a rather homogeneous gap func- 
tion in Sr 2 Ru04. Moreover, the action S'sd^i] in Eq- (§|) 
as a function of the charge, flavor, spin, and spin-flavor 
Bose fields is manifestly invariant under the mirror op- 
erations M. v [32] ]. The action including the interaction 
thus points towards a A\ g symmetry of the gap function. 
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This results from the fact that the charge, flavor, spin, 
and spin-flavor fields are linear combinations of the fields 
of the d zx and d yz orbitals. The A\ g symmetry is con- 
sistent with thermal conductivity measurements ^,[T~ 



well as with the geometry of the upper critical fields p4[ 
Point contact experiments also do not reveal any signif- 
icant in-plane anisotropy of the superconducting order 
parameter (|,|3j|. 

Finally, from the back-scattering terms in the action 
Ssg [4>i] the existence of two degenerate superconducting 
saddle points leading to two degenerate order parameter 
components can be deduced p3| |. Each component has 
a two-fold symmetry axis. Indeed, the existence of such 
two order parameter components with a slight spatial 
anisotropy in S^RuO is implied by the existence of two 
upper critical fields |24||39|. Comparison with the criti- 
cal field measurements suggests JIo| that the components 
are 93% isotropic. The two components are classified as 
flavor components f2f iX and Qf yV pjjfl . 

Since Clf tX and Qf tV are degenerate in the absence of 
fields breaking the ^-rotational symmetry the system 
can fluctuate between the two components in the ordered 
phase giving rise to a Goldstone mode jyj. This mode 
accounts for [^3| the gapless quasi two-dimensional exci- 
tations observed in the superconducting phase in various 
experiments P Jl5| , p^| . The presence of such a mode in 
the superconducting state finds support in the softening 
of the in-plane elastic constants recently observed in ul- 
trasonic measurements p2[ . 

In summary the notion of spin-triplet, orbital-singlet 
pairing in Sr 2 Ru04 leads to a straightforward physical 
picture that is consistent with fundamental experimental 
observations such as the sensivity to impurities, the sym- 
metry of the two upper critical fields and the fluctuations 
in the ordered phase. Unlike previous theories predicting 
p wave symmetry the two-component order parameter 
can be considered as an extended s wave with Ai g sym- 
metry and only slight anisotropy. The interactions play 
a crucial role in the in-plane correlations. Details of the 
non-perturbative approach and comparisons to experi- 
ments as well as the p-wave approach [35| are given in 
Refs. [11(230. 
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